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We describe the structural relations between nested harmonic sums emerging in the description of physical single 
scale quantities up to the 3-loop level in renormalizable gauge field theories. These are weight w=6 harmonic 
sums. We identify universal basic functions which allow to describe a large class of physical quantities and derive 
their complex analysis. For the 3-loop QCD Wilson coefficients 35 basic functions are required, whereas a subset 
of 15 describes the 3-loop anomalous dimensions. 



1. Introduction 

Scattering cross sections in renormalizable Quan- 
tum Field Theories which depend on a single kine- 
matic or mass scale obey a particularly simple 
form if dealt with in Mellin space. Contrary to 
the case in momentum-fraction space, where the 
variable z = p/P is referred to, or a representa- 
tion depending on z = s/t, with s and t Mandel- 
stam variables, the results are expressed in terms 
of nested harmonic sums with an outer summa- 
tion index N [1] . In this way a unique language is 
available in which the integrals emerging in higher 
order calculations can be expressed. Although it 
is not expected that even in this case the rep- 
resentation will be sufficient, cf. [2], it holds for 
all 3-loop calculations having been performed so 
far [3]. Furthermore, this representation holds in- 
dependent of the processes considered. 

The nested harmonic sums are recursively de- 
fined by 



N 



SbdN) 



j:'^^s,ik) 



k=l 



k\b\ 



(1) 



where a — (ai,...,a;) and b,ai G ^\{0} . The 
weight w and depth d of a harmonic sum Ss{N) 

are given by w = X]fe=i I'^fcl ^^'^ d = L At each 
weight w there are 2 • S™^""" harmonic sums. In 
the limit iV — > oo the nested harmonic sums turn 
into multiple zeta values Cs, [4]. As well known, 
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there is a large number of algebraic relations be- 
tween the multiple zeta values due to the shuffle- 
and stuffle relations, cf. [5]. The nested harmonic 
sums form a quasi-shuffle algebra, [6], with the 
shuffie product LU , which maps two index sets 
in all possible orders preserving the order of the 
original sets. The algebraic relations of the har- 
monic sums are well known, [7] , and are related to 
the shuffle product and the structure of the index 
sets. The number of the respective basis elements 
are counted by the Lyndon words or due to the 
Witt relations [8]. Let us consider the harmonic 
sums up to weight w=6. The reduction obtained 
cumulatively by the algebraic relations #r com- 
pared to the original number #c is 
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#c 
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26 


80 
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728 


#r 
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13 


31 


79 


195 



In the physical applications we observe [2, 9, 10] 
that one may find representations in which the 
index { — 1} never appears in a nested harmonic 
sum. The number of harmonic sums of this type 
iV_,|_i}(w) and the corresponding number of ba- 
sis elements N^'f^^-^y{w) are given by 
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with /i(d) the Mobius function. This leads to the 
following reduction 
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#c 
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11 


28 


69 


168 


#r 


1 


3 


7 


14 


30 


60 



Note that already the numbers at input are lower 
than the numbers after the algebraic reduction in 
the former case. 

All further relations between the nested har- 
monic sums are called structural relations. They 
result from the mathematical structure of these 
objects beyond that given by their indices. 

2. Structural Relations 

In the physical application nested harmonic sums 
emerge as a consequence of the light-cone expan- 
sion [11] and related techniques, observing the 
crossing relations for the respective processes, cf. 
e.g. [12]. The physical quantities are thus defined 
at either even or odd integers N. To use these 
expressions in experimental data analyzes it is, 
however, necessary to map these representations 
back to momentum-fraction space. This requires 
an analytic continuation of the nested harmonic 
sums from the even, resp. odd, values of N to 
N E C Eventually one has to derive the com- 
plex analysis for the nested harmonic sums. In 
passing, N takes values N G Q and N E M., 
which leads to new relations, as we will outline 
below. From a practical point of view, the most 
complicated part consists in deriving the corre- 
sponding representations for N E C for a large 
number of harmonic sums. Whenever possible, 
we will seek equivalence classes for these sums, 
the elements of which can be easily accessed ap- 
plying some operator, which can be straightfor- 
wardly realized even in the final numerical preci- 
sion representations. 

Harmonic sums can be represented in terms of 
Mellin integrals, cf. [1], 

SsiN) = M[fs{x)]iN) ^ f dx^f^\,,gx''-' . 

Jo ± ± X 

The first structural relation one obtains at w=l 
by decomposing 
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(3) 



which yields 
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f3{N) 
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(4) 



Since S'_i(iV) = (-l)^/3(iV-h 1) -ln(2), S'i(7V) = 
^{N + 1) -I- -fE, the sum S-i{N) is not alge- 
braically independent of Si{N) for N e Q. At 
higher weights ([3]) can be generalized using 



2k- 



ln(l-a;2) = 
-Lifc(a;2) = 



ln(l -x)+ ln(l + x) (5) 
hikix) +Liki~x) (6) 



as numerator function. 

Considering N E R one may define the differ- 
entiation of harmonic sums via 



— SsiN) = MMx) Mx)]{N) 



In this way we obtain 
d 



S2iN) 



dN 



Si{N) + C2 



(7) 



(8) 



In general one observes that the set of harmonic 
sums, extended by the multiple zeta values, is 
closed under differentiation. In particular one 
may form equivalence classes through differenti- 
ation and consider the harmonic sum with the 
lowest weight of this class as its representative. A 
consequence of these two relations is, that all sin- 
gle harmonic sums may be traced back to Si (N) . 
The complex analysis of the whole class derives 
from that of S'i(A^). 

A third class of structural relations emerges 
through the respective iterated integral represen- 
tation [13], combined with algebraic relations. At 
\N—2 one obtains 



M 



ln{l 



[N) = -M 



ln(l + x) 



1 + x \ L 1 + a; . 

mN)+jE+lnmi3iN)+P'{N) 



(N) 
(9) 



In this way we identify the first basic function 
ln(l + x) 



Fiix) 



(10) 



1 + x ' 

which is related to the harmonic sum S-i,i{N). 
Although in physics problems the index { — 1} is 



not emerging, Fi{x) is still useful to express a 
series of harmonic sums. All w=2 harmonic sums 
reduce algebraically or are related to Fi{x) and 
single harmonic sums. 

At w=3 we apply ([6]) for the first time. To- 
gether with a relation similar to ([9]) all sums 
which do not contain {—1} as index reduce to 
the Mellin transforms of 
Li2(a;) 



^2,3(2:) 



(11) 



l±x 

Let us now study the double sums in general. 
Using the above relations one may show that 
Nielsen integrals [14] are sufficient to express all 
double sums in terms of Mellin transformations 
of the functions 

Life (a;) 



FkM^) 



l±x 



(12) 



For even weight w there is an Euler relation and 
only Lifc(a;)/(1 + x) contributes. Let us illustrate 
this for 52,3 (iV). One derives 



S2AN) ^ KiSiiN) + 



M 



ln(a;) [S'i,2(l - xXs] + 3[Si,3(l - x) - (4] 



(N) 



where the functions 5*1.^(1 — x) are polynomials 
of ln(a;),ln(l — x) and Life(x). Therefore S2,3{N) 
is related to S'4_i(7V) up to derivatives of known 
harmonic sums of lower weight. 

In the following we only have to consider sums 
with depth d > 3. We derive all one-dimensional 
integral representations for the harmonic sums up 
to w=6 which do not contain an index {—1} and 
use the above relations. We obtain the following 
set of basic functions, [2,15]. 
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S,,3{x)/(x±l) 
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where 
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(14) 



and Na{x)\^^^ 3 are polynomials of Nielsen in- 
tegrals. Up to w=5 the numerator functions are 
Nielsen integrals, while at w=6 some of the nu- 
merator functions are general harmonic polylog- 
arithms over the alphabet {0,1,-1} [16] in the 
representation we derived. Yet it may still be, 
that an equivalent representation can be found 
over a two-letter alphabet, which will be investi- 
gated further. |3 

The Mellin transforms of the basic functions in 
(|13|) span the space of the harmonic sums which 
occur in higher order calculations for single scale 
quantities in Quantum Chromodynamics. Ex- 
amples are the anomalous dimensions and Wil- 
son coefficients in deeply-inelastic scattering up 
to 3-loops [3, 18]. With rising order in the cou- 
pling constant ag the number of contributing ba- 
sic functions is given by 

0{as) Wilson coeff./anom. dim.# 1 
0{ag) anomalous dimensions # 2 



^Structures of this type have been found at least in case 
of the multiple zeta values recently [17]. 



0{al) Wilson coefficients # < 5 
0{al) anomalous dimensions^ 15 
0{al) Wilson coefficients # 35 



(15) 



We represented a wide class of massless and 
massive 2-loop quantities in terms of harmonic 
sums, among them the unpolarized and polar- 
ized Drell-Yan cross section and the hard scatter- 
ing cross sections for hadronic Higgs-boson and 
pseudoscalar-boson production in the heavy top 
quark limit [19], the unpolarized and polarized 
time-like anomalous dimensions and Wilson co- 
efficients [20], the polarized anomalous dimen- 
sions and Wilson coefficients [21], the heavy fla- 
vor deep-inelastic Wilson coefficients in the limit 
Q^ ^ TO^ [22], including the 0{e) terms, as well 
as the virtual- and soft corrections to Bhabha 
scattering [10]. All these quantities fall into the 
third class above and are represented by the re- 
spective basic functions and simple polynomial 
factors in the Mellin variable. The same struc- 
ture is obtained for other similar processes more. 

3. Complex Analysis 

Finally we have to derive the analytic continu- 
ation of the nested harmonic sums iV S N — *■ 
iV G C. In the past precise numerical representa- 
tions were derived in [23] . |j Here the integrand is 
dealt with using the MINIMAX-mcthod [25], which 
leads to an adaptive representation in the com- 
plex plane in terms of a rational function. We 
seek, however, for a rigorous representation. 

It is well-known [26] that Mellin transforms of 
the type 



n{z) 



ifiii-t) 



dt t^^^ ip{t); 



^akt'' , 



k=0 



are factorial series, 



n{z) 



Ofc+ifc! 



^^ z(z+l). 



. (z -f fc) 



(16) 



(17) 



^This is also possible in the case of the heavy quark Wilson 
coefficients, including power corrections, cf. [24]. 



(f>(i) has to be analytic at t = 1. Vt{z) is a mero- 
morphic function for z € C, with poles at the 
non-positive integers. It obeys a recursion for 
z — > z -|- 1 and has an analytic asymptotic rep- 
resentation. The choice of the basic functions in 
(I13p is not immediately suited in all cases, as some 
of the functions possess branch points at t = 1. 
However, one may map this behaviour using re- 
lations for the basic functions, as for the change 
of the argument t ^ (1 — t). Let us consider the 
example 

F^iN) = M[Li2(a;)/(l -f x)\{N) . (18) 

The recursion relation is given by 



Fsiz + l) = -Fsiz) 



C2 



V'(z + 1) +7B 



(19) 



We map 



Li2(z) ^ -Li2(l -z)- ln(z) ln(l ~ z) + (2 , (20) 
and derive the asymptotic representation for 
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Li2(l-z) 
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1 
4iV3 


7 1 
24 7V4 


1 1 
3 7V5 


73 1 

^ 120 iV6 ■ ■ ■ 



(21) 

We stifl need to consider M[ln(z)ln(l - z)/(l + 
z)]{N), which is a derivative of the function 
Fi (N) and a few simpler terms. Representations 
like this can be obtained for all the basic func- 
tions. In this way analytic expressions for the 
Mellin transforms, resp. nested harmonic sums, 
which are needed up to the level w=6, are ob- 
tained. The above relations can be tuned to any 
numerical accuracy by analytic means. 

4. Conclusions 

In perturbative higher order calculations in renor- 
malizable Quantum Field Theories nested har- 
monic sums form an appropriate way of represen- 
tation for single-scale quantities, at least up to 3- 
loop orders. The number of multiple nested har- 
monic sums grows exponentially with the weight 



w. The harmonic sums obey algebraic and struc- 
tural relations, which reduce the number of ob- 
jects to a small set of basic functions. The quanti- 
ties are meromorphic functions in C with poles at 
the non-positive integers, obey recursion relations 
and possess an analytic asymptotic representa- 
tion. Due to this these functions are known in an- 
alytic form. The basic functions emerge as unique 
quantities in a large variety of higher order calcu- 
lations and hint to a rather simple mathematical 
structure behind Feynman diagrams. The sim- 
ple pattern, known in case of zero-scale quanti- 
ties, where only a few basis elements span the 
space of all multiple zeta values, is rediscovered 
for the next more complicated class of processes 
and found to be of quite similar structure. It is 
an observation that harmonic sums with an index 
{ — 1} can be avoided in the physical quantities, 
for yet unknown reason. Due to this the complex- 
ity of 728 objects at w=6, anticipated originally, 
reduces to 168. The latter ones can now be repre- 
sented by 35 basic functions, which is an essential 
compactification . 
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